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We report theoretical calculations of the spectrum of third sound modes for 
a cylindrically symmetric film of superfluid 3 He, and compare these results 
with experimental data for the mode frequencies and amplitude spectrum of 
surface waves of superfluid 3 He films. 
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1. INTRODUCTION 

Third sound is an oscillation of the superfluid component in the plane of 
the film, with the normal component clamped by viscosity to the substrate 
that supports the film. This sound mode was first discussed theoretically by 
Atkins for superfluid 4 He films in the context of the two-fluid modelP The 
oscillations of the superfluid component cause the level of the film to locally 
rise and fall (Fig. QJ. The restoring force that stabilizes the oscillations is 
the van der Waals attraction between the helium and the substrate. The 
possibility of observing third sound in superfluid 3 He films was discussed 
by Eggenkamp et al., 2 and the first experimental report of surface modes 
identified as third sound was by Schechter et al. 3 Here we report theoretical 
calculations of the spectrum of cavity modes of third sound in superfluid 3 He 
films with cylindrical symmetry and compare the results with the observed 
spectra. 



2. CAVITY MODES 



We consider 3 He films with thicknesses of order d ~ 300 — 400 nm, 
roughly 3 — 6 coherence lengths at zero pressure. These are thin films in 
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Fig. 1. Third sound oscillations of the superfluid produce surface waves. 

the sense that they are expected to be in the either the axial or planar 
phase. 4 In either case the orbital anisotropy axis, is anchored normal to 
the substrate. The analysis reported here assumes that the t vector remains 
fixed in the presence of small oscillations of the surface. We expect this to 
be a good approximation for long-wavelength, low-frequency third sound, 
since the coupling of the low-frequency normal-flapping mode in 3 He-A to 
transverse currents vanishes in the limit of exact particle-hole symmetry.^ 
We also neglect edge currents associated with the broken chiral symmetry 
of 3 He-A, which may play an important role in nearly two-dimensional films 
with specular surfaces, but are not relevant for third sound in the geome- 
try considered herePThus, in zero magnetic field the hydrodynamics of the 
superfluid film is defined by the phase mode and the surface excitation of 
the film, and reduces to scalar hydrodynamics. At temperatures T ~ mK 
the viscous penetration depth is of order S = \/2r]l pu « 1 cm ^> d, for low 
frequency third sound [y ~ 1 Hz) (77 and p are the viscosity and density of 
normal 3 He), and so we can safely consider the excitations (normal com- 
ponent) to be viscously clamped to the substrate. Furthermore, for third 
sound, the motion of the condensate is incompressible and confined to the 
plane of the film. Oscillations of the condensate then lead to height oscilla- 
tions of the film. The hydrodynamic equations governing the excitation of 
the film include the continuity equation for mass conservation, 

dC 

^ + ^V-v s = 0, (1) 

where p is the mass density, ( is the deviation of the film height from its 
equilibrium value, d, and p s is the superfluid density averaged over the thick- 
ness of the film. The pair momentum, or superfluid velocity, related to the 
phase 1} of the condensate's wave function, p s = |V$ = 7713 v s , describes 
the in-plane motion of the condensate and obeys the Josephson equation of 
motion, 

^ = -V M = -/VC + s/VT, (2) 
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where p is the chemical potential, / is the van der Waals force per unit 
mass, Sf is the entropy per unit mass and T is the absolute temperature. 
The oscillations of the condensate are coupled to entropy oscillations, and 
thus the height oscillations of third sound are coupled to the heat current in 
the film. This coupling, described by the heat transport equation, 

8T dC 

PCfd^ + ps f Tf t - dn f V 2 T = , (3) 

where Cf is the specific heat and Kf is the thermal conductivity per unit 
mass, is a key mechanism for damping of third sound cavity modes. 

If we neglect heat transport in the film, then the height deviation, £, 
and the superfluid velocity, v s , obey two-dimensional wave equations with 
the phase velocity given by c = yj p s (T)df/p. For cylindricalry symmetric 
films the eigenfunctions are products of Bessel functions (regular at r = 0) 
and azimuthal harmonics, J m (kr) e im ^ , where k = cj/c is the wavevector 
for the mode and m = 0, ±1,±2, .... The cavity mode eigenfrequencies 
are determined by the boundary conditions at the edge (r = R) of the 
film. For fixed boundary conditions, ((R) = 0, the resonant frequencies and 
wavevectors are u mn — c(T) a mn /R, where a mn is the n th zero of J m {x). 
For free boundary conditions, d r ((R) = 0, u mn — c{T) a f mn /R, where (a f mn ) 
are the zeroes of the derivative, J^x). The general solution to the wave 
equation can then be expressed as an eigenmode expansion, 

C(r, <f>, t) = ^ A mn J m (oj mn r/c) e im <t> e"*"™ * , (4) 

where the amplitudes, A mn , are determined by driving forces that couple to 
the displacement. 



3. ANALYSIS 

Schechter et alP observed standing waves of third sound on a circular 
substrate by exciting the surface with an a.c. electric field normal to the 
film. As noted in Ref. 7 the modes for thicker films (250 — 400 nm) are 
reasonably well indexed by assuming fixed boundary conditions, whereas 
thinner films (d < 250 nm) are better indexed assuming a free boundary 
condition. The mode frequencies for a 250 nm thick film are shown in Fig. [2j 
The temperature dependence reflects that of the superfluid density, which 
is approximately described by the Ginzburg-Landau form, p s (T)/p oc 1 — 
T/T^ lm . The fit of the theoretical mode frequencies to the experimental data 
is done assuming that only the lowest frequency modes with m = 0, 1 are 
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Fig. 2. Mode frequencies for a 250 nm film (symbols) from Ref. 7, The solid 
curves are the calculated mode frequencies with mode indices as labelled for 
fixed boundary conditions. 



detected. In general there are higher frequency modes that may contribute, 
as well as modes with m > 2; these modes are also shown in Fig. [21 for 
comparison. It is not apriori obvious that these modes can be ignored. 
There are significant differences between theory and experiment, particularly 
for the lowest order modes. In order to more precisely identify the modes 
and to test the hydrodynamic theory for third sound oscillations of 3 He 
we have calculated the amplitude response of the film. The calculation 
includes the sensitivity of the detector for the various modes, and provides 
a determination of the amount of energy that is pumped into each mode. 

The excitation field, or driving force, enters via the electrostatic energy 
density of the thin film which is a weak dielectric. The electro-chemical 
potential is then 

E 2 e — 1 

p = no(T,p) - , (5) 

07T p 

where p is the mass density, e is the dielectric constant of the liquid 3 He 
and E is the electric field. The response of the film to an a.c. voltage, 
V(t) = Vq + V\e lujt , is governed by the wave equation with the driving force 
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on the right-hand side, 

- C 2 V 2 C = -c\T)A Y,fmnJ m (a mn ^)cos(m<f>) e*"* (6) 

where A ~ 2A is of the order of the surface wave amplitude and f mn is the 
strength of the excitation force for mode (m,n). The capacitive excitation 
ring is displaced above the film and is slightly tilted about an axis (y-axis) 
the horizontal plane of the film. 3 If we measure the azimuthal angle with 
respect to the orthogonal direction in the plane (x-axis) then the resulting 
forces on the film are symmetric with respect to <j) —> — <j). The surface 
displacement is then 

C(r, <f>,t) = A c 2 V J mn 2 J m (a mn r/R) cos(m^) e*"* . (7) 



The capacitive detection method used in Refs. EE allows one to measure 
the average displacement of the film surface. The amplitude spectrum, A(uj), 
is proportional to the average displacement and is given by 



A{ui) =A[ 2 -^) c\T) 



F 

x mn 



u mn ~ u2 
m.n 111,1 



(8) 



where F mn is the product of the excitation force (f mn ) and the detector 
sensitivity (D mn ) for mode (m,n): F mn — f mn x D mn . The detector sen- 
sitivity is determined by the electric field distribution for the geometry of 
the capacitor, which we calculate for the experimental setup of Ref. El The 
modes shown in Fig.[3|for the response are indeed the m — 0, 1 modes, which 
justifies the earlier assumption. However, the spectra shown in Fig. [HI differ 
significantly from the experimental spectra. There are many more modes 
visible in the calculated spectrum. This is a consequence of the assumption 
of frequency independent damping for the modes; constant damping does 
not suppress the higher frequency modes. 

In general the damping of third sound is frequency dependent. In the 
following we calculate the frequency dependent damping rate for third sound 
that results from heat transport induced by the oscillations of the entropy 
associated with oscillations of the condensate. The resulting amplitude spec- 
trum is obtained by replacing u — » uo + ir(T, uj) in Eq. [HI in the response 
formula. The damping rate contains two terms, V — T$ + Ti(T, cj), where Tq 
is a constant damping rate which we attribute to roughness of the substrate 
and Ti(lj, T) is the damping obtained from the heat transport equation (Eq. 



2 \c f J fdpc 2 {T) l-t 
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Fig. 3. Spectral weight for various modes of third sound folded with the 
detector sensitivity. The spectra are calculated assuming a constant, fre- 
quency independent damping, T, with free boundary conditions (left) and 
fixed boundary conditions (right). 



where t = T/T^ lm . Note that the damping rate diverges for T — > T^ lm be- 
cause the superfluid stiffness vanishes whereas the thermal conductivity is 
finite. We have also assumed that surface pair breaking yields superfluidity 
with a gapless spectrum, in which case Sf = Cf « 7 S T and Kf oc TP The 
resulting damping rate leads to strong suppression of the higher frequency 
modes and much closer agreement with the experimental spectra. The cal- 
culated results for the amplitude spectra with both fixed and free boundary 
conditions are shown in comparison with the experimental spectra in Figs. 
0] and [3 For the thick film with d — 380 nm the theoretically calculated 
spectrum are in reasonable agreement with the experimental spectrum for 
temperatures T < 0.64 mK. At higher temperatures mixing or crossing of the 
modes is seen in the experimental spectrum. Such features are not observed 
in the theoretical spectra obtained for either type of boundary condition. 
The discrepancies are more severe for thinner films. In this case neither 
boundary condition yields agreement between the calculated and observed 
spectra. Experimentally, there is a much richer mode structure than the 
theoretically predicted spectrum. 



4. CONCLUSIONS 



In summary, the thicker films (250-400 nm) at low temperatures are 
reasonably well described by the equations of two-fluid hydrodynamics with 
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Fig. 4. Experimental spectra from Ref. 0for a "thick" film with d = 380 nm 
(left) and the corresponding theoretical spectra calculated for fixed bound- 
ary conditions. Spectra for different temperatures are shifted vertically as 
indicated in the experimental figure. 




Fig. 5. Experimental spectra from Ref. 0for a "thin" film (230 nm) and the 
corresponding theoretical spectra calculated with free boundary conditions. 
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damping due to surface roughness and thermal transport. The latter mech- 
anism suppresses the higher frequency modes and diverges for T — ► T c film . In 
this regime the theoretically predicted modes and their spectral weights cor- 
respond to those observed in the experiment. The apparent mode splitting 
or mixing that onsets at higher temperatures, T > O.GST^ 1 ™, may indicate a 
cross-over to a another mechanism for dissipation or, what seems more likely 
given the abruptness of the onset, a phase transition of the superfluid phase 
of the film. For thinner films (170-250 nm) we are not able to describe the 
observed spectra using two-fluid hydrodynamics with either set of boundary 
conditions and the same damping mechanisms. Although the frequencies 
of the modes can be associated with modes calculated with free boundary 
conditions, the spectral weights for the modes calculated from the theory do 
not agree with experiment. This fact, and mixing or splitting observed in 
thicker films may indicate that there is a phase transition at temperatures 
well below T c for films with d <~ 3 — 6£o- 
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